Nonviscous motion of a slow particle in the dust crystal under microgravity conditions 
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Subsonic motion of a large particle moving through the bulk of a dust crystal formed by negatively 
charged small particles is investigated using the PK-3 Plus laboratory onboard the International 
Space Station. Tracing the particle trajectories show that the large particle moves almost freely 
through the bulk of plasma crystal, while dust particles move along characteristic a-shaped pathways 
near the large particle. In the hydrodynamic approximation, we develop a theory of nonviscous dust 
particles motion about a large particle and calculate particle trajectories. A good agreement with 
experiment validates our approach. 
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I. INTRODUCTION 

Complex or dusty plasmas [1, 2] are low temperature 
plasmas consisting of electrons, ions, neutral gas atoms, 
and charged microparticles or dust particles. Due to the 
higher electron mobility, these particles acquire large neg- 
ative charges up to 10 4 e, where e is the electron charge. 
Thus, the particle ensemble is a strongly coupled classical 
Coulomb system, in which plasma crystal, i.e., ordered 
structures of dust particles can be observed [3-8]. In ex- 
perimental setups, dusty plasmas are usually studied in 
gas discharges at low pressures, e.g., in radio frequency 
(RF) discharges [3] or in striations of dc discharges [9-11]. 
In ground-based experiments, the electric field present in 
a discharge compensates the force of gravity and plays 
the role of a trap for dust particle structures. Three- 
dimensional structures obtained in RF and dc discharges 
are inhomogcneous ones. A large homogeneous bulk of 
dusty plasma, which almost fills the entire discharge vol- 
ume, can be observed under microgravity conditions ei- 
ther in parabolic flights [12-16] or onboard the Interna- 
tional Space Station (ISS) [12, 17-21]. 

Recently, self-excited dust-density waves caused by 
streaming ions were studied [15], and freezing and melt- 
ing of dusty plasma upon variation of the neutral gas 
pressure was investigated [19]. The motion of a large 
particle inside the dust crystal of smaller particles was 
reported in [13, 16, 21]. It was demonstrated that su- 
personic motion is accompanied by the formation of the 
Mach cones, causing noticeable perturbation of plasma 
crystals. There were also relatively slow particles ob- 
served, moving with almost a constant velocity or even 
accelerating, which is incompatible with common notion 
of the interaction between the particle and dust crystal. 
Such particles appear to move freely, without long-range 
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distortion of the ambient structure. 

Using the PK-3 Plus laboratory onboard the ISS we 
have investigated motion of the dust particles in the 
vicinity of moving large particle. Such relatively slow 
particles are accelerated spontaneously on the outside of 
the particle cloud in RF discharge, penetrate into the 
dust cloud, and move inside it. We have made a series 
of snapshots of this motion in the plane of a laser sheet 
with a high-resolution camera, which allowed us to trace 
the trajectories of individual particles. Motion of small 
particles caused by the passage of a large particle has 
already been studied and used as a source of valuable 
information in the course of a ground-based experiment 
[22], where heavy particles were dropped down to the 
cloud of dust particles levitating in the plasma sheath re- 
gion of a RF capacitive planar discharge. Trajectories of 
particles suspended in low pressure glow discharges have 
also been studied in [23] . Dust particles were arranged in 
chains, and a heavy particle falling between neighboring 
dust chain bundles caused the elliptical motions of the 
background dust particles. 

Superposition of successive frames recorded in our ex- 
periment reveals that many more dust particles in a close 
neighborhood of a moving particle circumscribe typical 
a-shaped trajectories, while other particles remain al- 
most at rest. We interpret this motion as a nonviscous 
flow about a large particle, for which we used a hydro- 
dynamic approach. This approach is widely used as ap- 
plied to the collective motion of dust particles [2, 24-26]. 
Based on classical solution for the velocity field [27], we 
integrate it to obtain streamline pathways for individ- 
ual dust particles. A good agreement between recorded 
and calculated trajectories validates our basic assump- 
tion. Since a nonviscous motion is associated with zero 
momentum transfer, the force acting on a large parti- 
cle vanishes in the idealized hydrodynamic model. If the 
particle interact only with the dust particles in its close 
neighborhood, as in the case of a nonviscous flow, the 
estimate of the direct particle-dust interaction force sug- 
gests that it is not greater than other major forces oper- 
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ating in low-pressure gas discharges (such as the ion drag 
and electric forces). 

The paper is organized as follows. In Sec. II, we de- 
scribe the experimental setup and results of tracing the 
dust particle motion caused by the passage of a large par- 
ticle. A hydrodynamic approach is developed in Sec. Ill, 
where both numerical and analytical solutions for dust 
particle trajectories are obtained. Results are summa- 
rized in Sec. IV. 



II. EXPERIMENTAL SETUP AND RESULTS 

Here, we present an observation of a subsonic nonvis- 
cous motion of a large particle in the bulk of a three- 
dimensional complex plasma in the PK-3 Plus labora- 
tory onboard the ISS. Details on the setup can be found 
in [20] . The heart of this laboratory consists of a capac- 
itively coupled plasma chamber with circular electrodes 
of 6 cm diameter and 3 cm apart (see Fig. 1). Microparti- 
cles can be injected into the main plasma with dispensers. 
The particles then form a cloud around the center of 
the chamber, typically with a central void caused by the 
ions streaming outwards. Generally, some larger particles 
are present in the chamber as well. The origin of these 
particles is not yet understood; these might be, for in- 
stance, agglomerates or larger particles left over from pre- 
vious experiments (not removed during the cleaning pro- 
cedure). They normally accumulate themselves at the pe- 
riphery of the particle cloud, because of the dependence 
of the ratio of the electric and ion drag force on the par- 
ticle diameter [28]. Sometimes, one of these larger par- 
ticles gets sporadically accelerated and penetrates into 
the cloud - we shall term it projectile. The reason for 
such behavior remains a puzzle (for instance, it might be 
caused by a laser- induced rocket force [29]). Neverthe- 
less, these particles can be used as natural probes of the 
microparticle cloud [30, 31], in the same way as parti- 
cles injected on purpose [13, 16, 32]. Projectiles enable 
various active experiments to be performed in a wide 
range of testing parameters, thus providing the opportu- 
nity to reveal new features of strongly coupled complex 
plasmas. In particular, velocities of projectiles can vary 
significantly, and when they move through a cloud of 
background microparticles at supersonic velocities they 
excite Mach cones [21]. On the other hand, obvious dis- 
advantages of experiments with projectiles is that their 
velocities cannot be controlled. 

In what follows, we discuss in detail one example of a 
slow (subsonic) motion of a projectile. The experiment 
was performed during the 13th mission of PK-3 Plus on 
the ISS. Argon was used as buffer gas at a pressure of 
10 Pa, and the main microparticle cloud was composed 
of melamine-formaldehyde particles with 2.55 fj,m diam- 
eter. The projectile was accelerated from the side and 
penetrated into the main cloud, moving almost horizon- 
tally from the left to the right (see Fig. 2a, phase I). To 
determine the size of the projectile, this experiment was 



compared with others performed with larger particles, 
which allowed us to conclude that the projectile radius 
was a p — 7.5 /im. During its path towards the void, 
the velocity of the projectile varied, decreasing from 80 
mm/s to 37 mm/s. Inside the void, the horizontal mo- 
tion further slowed down. The particle was then accel- 
erated upwards (see Fig. 2a, phase II), where it again 
penetrated into the microparticle cloud (see the discus- 
sion of microparticle trajectories in Ref. [21]). In the re- 
gion above the void, the particle was slower than before. 
It did, however, push the microparticles away to clear its 
path. Figure 2b shows the projectile motion during phase 
II. The microparticles that were pushed away moved in 
vortices around the probe particle. The projectile accel- 
erated while traveling upwards through the microparticle 
cloud, so that its velocity increased from 7 to 14 mm/s. 
This value is still lower than the speed of sound. 

In the experiments, we monitored the motion of dust 
particles in the vicinity of the projectile, by manually de- 
termining its position recorded by a high-resolution cam- 
era at 50 frames/s (a video frame is illustrated in Fig. 2b). 
The horizontal and vertical resolution of the camera was 
11.3 /im per pixel and 10.3 /urn per pixel, respectively 

Experimental results on the motion of dust particles in 
the neighborhood of a slowly moving projectile are sum- 
marized in Fig. 3. The projectile motion from the void 
to chamber wall is represented by frames from bottom 
to top in Fig. 3a. Dust particles form a dust crystal, 
which can be represented as a set of the Wigner-Seitz 
cells around each particle. We estimate the dust par- 
ticle number density as rid = 3 x 10 5 cm~ 3 , hence the 
cell radius is f d = (3/47rn £i ) 1 / 3 = 9.3 x 10~ 3 cm. We 
can define the coupling parameter of interaction between 
the projectile and dust particles (or the scattering pa- 
rameter) (3dp = 2Z p Zde 2 J 'XMdU 2 as the ratio of charac- 
teristic Coulomb energy Z p Zde 2 /X at the plasma (ion) 
screening length A (where Z p and Z d are the charges of 
the projectile and dust particle, respectively, in units of 
the electron charge e) to the kinetic energy of a dust 
particle M d u 2 /2 in the reference frame co-moving with 
the projectile (where M d = 1.31 x 10" 11 g is the dust 
particle mass and u ~ 1 cm/s is the projectile veloc- 
ity). Based on the analysis of recent melting experiments 
performed under identical conditions [19], we estimate 
A ~ 6 x 10~ 3 cm, Zd — —1200 and, assuming a linear 
scaling of charge with size, Z p ~ —7000. This yields 
Pdp — 100 and, hence, the estimate for the radius of 
strong interaction between the projectile and dust par- 
ticles gives [8] ~ A In ft dp — 2.7 x 10~ 2 cm. This value 
turns out to be fairly close to the average distance be- 
tween the projectile and nearest dust particle neighbors, 
R ~ 3.3 x 10~ 2 cm (see Fig. 3a). Figure 3b shows suc- 
cessive phases of motion of an individual dust particle. 
Trajectories of such particles are revealed by frame super- 
position (Fig. 3c). It is evident that many more particles 
that interact with the projectile move along loops of sim- 
ilar form. The absence of loops in some regions can be 
accounted for by the following facts: first, the projectile 
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trajectory is not exactly parallel to the laser illumination 
plane as well as dust crystal planes, so one can see only 
projections of dust particle trajectories. Second, visible 
pathways of the dust particles terminate at the points 
where particles go out of the laser sheet. Third, as the 
distance from the projectile increases, dust particles come 
at rest and become a part of the crystal. Upon such local 
crystallization, the particle can be pushed away from its 
initial trajectory. 

III. THEORY OF NONVISCOUS DUST 
PARTICLE MOTION 

Let us first estimate the major forces acting on the 
projectile, in order to enlighten the physics of interac- 
tion between the projectile and the surrounding complex 
plasmas. We shall confine ourselves to treatment of the 
neutral drag force F n , the ion drag force Fi, the electric 
force F e , and the force from dust particles Fd- 

We start the analysis with the neutral drag force, for 
which we can provide the most accurate estimate. This 
force is naturally pointed in the direction opposite to the 
projectile velocity. Assuming, as usual, diffuse scattering 
of neutrals from the particle surface (with full accommo- 
dation) we obtain [8] 

8\/ 27T o _c 

F„ = S — - — m n n n v Tn a p u ~ 3 x 10 dyn, 

where S ~ 1.4 is the accommodation coefficient, n n = 
2.42 x 10 15 cm" 3 is the number density of neutral gas at 
the pressure of 10 Pa, m n = 6.63 x 10 -23 g is the mass 
of argon atom, and vx n = iT n /m n ) x ^ = 2.50 x 10 4 cm/s 
is the thermal velocity of neutrals. 

To estimate the plasma-induced forces, we need to 
know the magnitude of a dc electric field, E. The field 
is generated due to ambipolar diffusion and can be de- 
duced from the scaling [33] eEL ~ T e , where L = 3 cm 
is the distance between the RF electrodes (see Fig. 1) 
and T e = 3.5 eV is the electron temperature (here and in 
what follows, the temperature is expressed in the energy 
units). Then the electric force is 

F e = \Z p \ eE = \Z p \ TJL ~ lO" 8 dyn. 

We note, however, that in dense clouds of dust particles, 
where the Havnes parameter is large, E can be dramat- 
ically reduced (by an order of magnitude or even more, 
see Refs [34, 35]). Therefore, the obtained value only be 
treated as a crude upper-bound estimate. 

The ion drag force is determined by the ion drift ve- 
locity Uj. Assuming the mobility-limited ion drift we get 
m = eE/mm ~ T e /(miVTin n o-i n L) ~ 3 x 10 4 cm/s, 
where i/j = VT%n n o~in is the ion-neutral collision fre- 
quency, VTi(= VTn) is the ion thermal velocity, and 
(Ti n — 10 -14 cm~ 2 is the ion-neutral collision cross sec- 
tion. The magnitude of the ion drag force is determined 
by the ion scattering parameter j3 p i = \Z p \e 2 /XTi, which 



is of the order of 10 for our conditions. In this case, 
the force on the projectile is dominated by the direct ion 
collisions (absorption) and can be estimated from the fol- 
lowing formula [8]: 

Fi ~ 7rA 2 In 2 ppimiUiuf ~ 10~ 8 dyn. 

Here, we took into account that the ion drift velocity is 
of the order of the thermal velocity, and estimated the 
ion number density rij from the quasineutrality condition 
at large Havnes parameter, n,i = \Za\ nd ~ 3 x 10 s cm~ 3 . 
Thus, F e and Fi occur to be of the same order of magni- 
tude. Note that the ion drag force, which always opposes 
the electric force, is pointed away from the center. 

As regards the interaction between the projectile and 
dust particles, we assume that the Coulomb interaction is 
screened at short distances. Then the upper-bound esti- 
mate would suggest total transfer of the momentum from 
dust particles that find themselves inside this cylinder to 
the projectile. The force is given by [8] 

F d = nR 2 n d M d u 2 ~ 10" 8 dyn, 

which is not greater than the upper-bound estimates for 
the plasma-induced forces. 

We conclude that the major known force acting on the 
projectile is the neutral drag which opposes the unknown 
driving force and, presumably, results in the observed 
(locally) steady-state motion (the neutral damping rate 
is F n /M p u ~ 10 s _1 , i.e., the transition to steady state 
occurs within ~ 0.1 s). Leaving the discussion of possible 
driving mechanisms aside, let us focus on the flow of dust 
particles caused by the projectile. We assume that the 
plasma crystal is melted in the neighborhood of the pro- 
jectile that is streamlined by dust particles. If the stream 
flow is a potential nonviscous one the total momentum 
is conserved, and Fd must vanish. This assumption can 
be validated by comparison between dust particle trajec- 
tories calculated theoretically and those obtained from 
experiment. 

We apply a hydrodynamic approach to describe the 
dust particle motion, which is widely used in physics of 
complex and dusty plasmas [2]. Based on the discus- 
sion above, we can treat the projectile-dust interaction 
as that of hard spheres characterized by the length scale 
R. Consider a steady flow of an incompressible liquid of 
dust particles, described by the Navicr Stokes equation, 

(v V)v = -V- +7yV 2 v, (1) 
P 

where v(r) is the velocity field, p and p are the pres- 
sure and density, respectively, and r\ is the shear viscosity 
of the particle fluid. For an incompressible fluid, when 
p = const, Eq. (1) is completed by the continuity equa- 
tion 

divv = 0. (2) 

Let us estimate the relevant Reynolds number, Re = 
puR/r/, which characterizes the relative importance of 
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the viscous term in Eq. (1). Using results of the 
molecular dynamics simulations of Yukawa fluids [36], 
we obtain that for our conditions the kinematic viscos- 
ity can be estimated as rj/p = (0.1 — 0.3)w p df^, where 
ujpd = ^AnZ'^e 2 nd/Md ~ 300 s _1 is the dust plasma fre- 
quency, so that the Reynolds number is Re ^10. Thus, 
we can reasonably neglect the viscosity and thus reduce 
Eq. (1) to the Euler equation for an ideal fluid, 



p(vV)v = — Vp. 



(3) 



For an irrotational flow of incompressible liquid (rot v = 
0), the substitution v = V<£> transforms (2) to the form 



0. 



(4) 



where ip is the velocity field potential. 

Consider a streamline about a sphere with the radius 
R moving with a constant velocity u. The boundary con- 
ditions arc v(oo) = and v(-R) • n = u • n, where r is the 
radius vector from the center of the sphere and n = r/r 
is the radial unit vector. The unique solution of Eq. (4) 
with these boundary conditions is ip(r) = — (R 3 /2r 2 )u-n, 
so that [27] 



v(r) 



2r 



3 [3n(u 



(5) 



The pressure field p(r) can be derived from Eq. (3) but 
it is unnecessary for our purposes. 

The streamlines can be obtained by integration of 
Eq. (5) in the laboratory frame of reference. Let the X- 
axis be directed along the motion of a sphere. Then vec- 
tor n has the components {(x — ut)/r, y/r, z/r}, where 

1/2 

r = Ux — ut) 2 + y 2 + z 2 ] and u • n = (x — ut)u/r. We 
introduce the variables C — (x ~ ut)/R, r) = y/R, and 
t = vt, where v = 3u/2R, to rewrite Eq. (5) as a set of 
differential equations, 



c- 
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£2 _|_ jy2 



d V , ,,2 



V 



T 

2N-5/2 



2N-3/2 
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3' 



(0) 



with the initial conditions £(0) = and r)(0) = 770 (|?yo| > 
1). Due to the axial symmetry, the third equation for 
Z -axis coincides with the second equation (6) and is 
therefore redundant. 

Because the decay of particle velocity (5) with the dis- 
tance from projectile is rather fast, |v| cx r~ 3 , we can 
assume that particles move in a thin fluid tube in the 
vicinity of a sphere with the radius R so that the differ- 
ence between r and R can be neglected and £ 2 + rj 2 
Equations (6) arc then reduced to 



1. 



dC 

dr. 



drj 



07, 



(7) 



with the initial conditions £(0) = and 77(0) = 1. Within 
the framework of this approximation, we have to assume 
also that a particle is quiescent until it finds itself on the 
surface of a moving sphere (projectile cell). This instant 
corresponds to the time t = —At and to the distance 
d between the particle and direction of projectile mo- 
tion (impact parameter). Obviously, solution of Eqs. (6) 
obeys the condition ^(±00) | = d/R. Due to the time 
reversibility of Eqs. (7), the particle must stop at t = At 
and further remain quiescent, the total time of motion be- 
ing 2At. Thus, the initial conditions should be completed 
by the 'collision' condition rj(At) = ±d/R, where the sign 
of d defines the direction of particle motion. We integrate 
Eqs. (7) to derive ((t) = — tanhr and r/(r) = (coshr) -1 . 
In conventional units, we obtain 



x 
R 



-vt — tanh vt, 
' 1 



'J 



R cosh vt ' 



(8) 



From the second equation (8) we obtain 



At = 




(9) 



If 1 — (d/R) <C 1, we find approximately At ~ 

1/2 — 

(2R/d— 2) for a grazing collision. At (d/R) <C 1, we 
have At ~ ln(2^/c?). Streamline (8) is shown in Fig. 4 
for d/R = 0.163, which corresponds to vAt = 2.5 in the 
framework of used approximation. It is seen in Fig. 4 
that the approximate solution almost coincides with the 
numerical solution of equation set (6) corresponding to 
the same d/R. Our approximation is valid even at a large 
impact parameter d = R because in this case, numerical 
solution of (6) yields y(0) = 1.33-R, i.e., the distance be- 
tween the particle and projectile is still not much different 
from R. Figure 5 shows a good correspondence between 
solution (8), (9) and experimentally observed traces of 
individual particles extracted from Fig. 3c in their closed 
portions. Indication of this correspondence is a good 
agreement between the ratio of the height of a closed 
loop to R determined experimentally (Figs. 3c and 5), 
which amounts to about 0.28, while the theoretical value 
is 0.277 (cf. Fig. 4). This ratio is independent of the loop 
position relative to the projectile path. Note that since 
the streamline takes place in a close neighborhood of the 
projectile and farther regions are crystallized, an approx- 
imate solution (8), (9) utilizing the same approximation 
seems to be even more adequate to treated system than 
the numerical solution. Open portions of particle traces 
reveal some differences from theory. In Fig. 5, the upper 
trace portion of the bottom loop on the left-hand side 
seems to disappear due to the fact that the particle go 
out of the illumination plane. On the contrary, the upper 
trace portion of the top loop on the left-hand side would 
not appear at all but the particle was probably pushed 
towards the projectile due to spatial re-distribution of 
particles in the course of local crystallization. 
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We have demonstrated that the stream flow of dust 
particles is almost a potential nonviscous one. In an ideal 
model, the momentum transferred from dust particles to 
projectile and back must coincide, so that the force Fd of 
interaction between the dust particles and projectile van- 
ishes. In reality, however, F d is never equal to zero. It fol- 
lows from solutions for the particle trajectories obtained 
above that the motion of particles outside the cylinder 
(whose axis coincides with the projectile trajectory) can 
be disregarded, as we originally assumed. The lower- 
bound estimate stems from the momentum not trans- 
ferred back to the projectile, due to collisions between 
neutrals and dust particles moving along the loops with 
the velocity estimated as u/2: 

where ad = 1.275 x 10~ 4 cm is the dust particle radius. 
Thus, Fd proves to be of the same order of magnitude 
as its upper-bound estimate obtained above and as F n . 
However, if treated damping effect was taken into account 
in Eq. (1) the solution for trajectories (8), (9) would not 
change because it is a consequence of unchanged Eq. (2) 
with the same boundary conditions. 



IV. CONCLUSION 

In this study, we have observed that a heavy particle 
(projectile) going through the dust crystal with a rel- 
atively low (subsonic) velocity appears to move almost 
freely. We have proposed a new explanation of this fact 
implying a nonviscous flow of the dust particles about 
the projectile, in which the force of interaction between 
the projectile and dust particles vanishes. To prove this, 
we monitored the particle dynamics by means of a high- 
resolution camera. Superposition of successive frames re- 
vealed typical a-shaped pathways of dust particles in the 
neighborhood of projectile. Since the system is strongly 
coupled, it can be divided into the Wigner-Seitz cells 
around each particle including the projectile, whose mo- 
tion does not break close ordering. However, the projec- 
tile melts the crystal in its neighborhood, and the dust 
particles flow about the projectile cell boundary in such 
a way that the distance R between individual dust par- 
ticles and the projectile, which can be measured using 
recorded frames, is the sum of radii of the Wigner-Seitz 
cells around dust particles and the projectile. In the 
framework of our model, we have to imply that the cells 
around particles rather than the particles themselves take 
part in the interaction. In the absence of direct interpar- 
ticle contacts, the shear viscosity must vanish and we can 
consider the nonviscous flow of the cells. 

The classical solution for the velocity field in a po- 
tential nonviscous flow about a sphere moving with a 
constant velocity can be regarded as a set of differential 
equations for the trajectories of individual dust particles. 



We have solved these equations numerically and obtained 
an analytical solution, which approximates closely the 
numerical one. This solution, which defines the coordi- 
nates of dust particles as explicit functions of time, is in 
a good agreement with particle trajectories obtained ex- 
perimentally by superposition of successive frames. This 
supports our main assumption concerning the nonviscous 
nature of particle flow. Because such motion implies no 
momentum transfer from liquid to the projectile, the drag 
force Fd arising from the interaction between the projec- 
tile and dust particles must vanish. Note that examples 
of such type of dust particle collective motion were pre- 
viously not encountered in physics of complex and dusty 
plasmas. 

A real system does not show idealized behavior. Not all 
the particles move along the a-shaped pathways, which 
we attribute (i) to the fact that the illumination plane 
does not coincides with one of the crystalline planes and 
is not exactly parallel to the direction of the projectile 
motion, (ii) to the finite thickness of the laser sheet, so 
that particles can leave the illuminated space, and (iii) 
to a local crystallization after the projectile passage. In 
addition, the force of interaction between the projectile 
and dust particles does not vanish exactly. However, a 
general conclusion of the developed theory, that only the 
dust particles inside a cylinder with the radius R are 
involved in collective motion, is strongly confirmed by 
Fig. 3a. Note that this is a consequence of the boundary 
condition imposed on the projectile cell, which implies a 
nonzero tangential velocity field component. In the case 
of a viscous flow (e.g., Refs [24, 26]) the situation can be 
different due to vanishing of this component, resulting to 
much larger number of moving particles. 

We can conclude that a large particle can move almost 
freely inside the bulk of a strongly coupled dust crystal. 
Our investigation also points to the fact that a hydrody- 
namic approach can be valid at small length scales down 
to several intcrparticlc distances. The processes of local 
dust crystal melting and freezing, which were not treated 
in detail in the present study, will be addressed in future 
work. 
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FIGURE CAPTIONS 

Fig. 1: The cross-sectional schematic of the plasma 
chamber 

Fig. 2: (a), snapshot of the microparticle cloud from 
quadrant view camera (field of view is 35.7 mm x 26.0 
mm) . On the left-hand side, a small track of initial move- 
ment of the probe particle (projectile) into the cloud is 
visible. The horizontal dashed line arrow shows schemat- 
ically the path of projectile on the left-hand side (phase 
I); the vertical dashed line arrow, the path of projectile 
from the void into the upper part of microparticle cloud 
(phase II). (b), snapshot of the microparticle cloud from 
high-resolution view camera (field of view is 8.1 mm x 
5.9 mm). A small track on the right-hand side shows the 
projectile motion in phase II 

Fig. 3: Recorded positions of dust particles and the 



projectile (negative images). Time interval between in- 
dividual frames is (a), 0.1 s and (b), 0.02 s; (c) is a super- 
position of 12 frames with the time interval of 0.02 s. The 
projectile motion from void to chamber wall corresponds 
to frames from bottom to top 

Fig. 4: Calculated trajectory of a dust particle mov- 
ing about the projectile for v/S.t = 2.5. Solid line shows 
calculation using (8), (9) and dashed line, numerical inte- 
gration of Eqs. (6) for the same d = 0.163-R. Dots mark 
the intervals vt = \/A 

Fig. 5: Individual dust particles moving along the 
a-shaped pathways. Traces are the superposition of 
frames obtained in the experiment (negative image, verti- 
cal trace marks the projectile trajectory), solid lines show 
calculation by formulas (8) and (9), dashed line indicates 
the trajectory portion 'forced' by local crystallization 
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